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ABSTRACT 

In this paper we continue our study of CMB TE cross correlation as a source of 
information about primordial gravitational waves. In an accompanying paper, we con- 
sidered the zero multipole method. In this paper we use Wiener filtering of the CMB 
TE data to remove the density perturbation contribution to the TE power spectrum. 
In principle this leaves only the contribution of PGWs. We examine two toy experi- 
ments (one ideal and one more realistic), to see how well they constrain PGWs using 
the TE power spectrum. We consider three tests applied to a combination of observa- 
tional data and data sets generated by Monte Carlo simulations: (1) Signal-to-Noise 
test, (2) sign test, and (3) Wilcoxon rank sum test. We compare these tests with each 
other and with the zero multipole method. Finally, we compare the signal-to-noise 
ratio of TE correlation measurements first with corresponding signal-to-noise ratios 
for BB ground based measurements and later with current and future TE correlation 
space measurements. We found that an ideal TE correlation experiment limited only 
by cosmic variance can detect PGWs with a tensor-to-scalar ratio r = 0.3 at 98% 
confidence level with the S/N test, 93% confidence level with the sign test, and 80% 
confidence level for the Wilcoxon rank sum test. We also compare all results with cor- 
responding results obtained using the zero multipole method. Wc demonstrate that to 
measure PGWs by their contribution to the TE cross correlation power spectrum in a 
realistic ground based experiment when real instrumental noise is taken into account, 
the tensor-to-scalar ratio, r, must be approximately four times larger. In the sense to 
detect PGWs, the zero multipole method is the best, next best is the S/N test, then 
the sign test, and the worst is the Wilcoxon rank sum test. 

Key words: cosmic microwave background - polarization - gravitational waves - 
cosmological parameters 



1 INTRODUCTION 

Primordial gravitational waves (PGWs) (tensor) generate 
negative temperature-polarization (TE) correlation for low 
multipoles, while primordial density (scalar) perturbations 
generate positive TE c orrelation for low mu lt ipoles (s ee 
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test based on this sig nature (the zero multipole method, see 
(2007)) is useful as an insurance against false 
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detection or as a monitor of imperfectly subtracted system- 
atic effects. 
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In this paper, we analyze an alternative method. This 
method uses Wiener filtering to remove the contribution 
of density perturbations to the TE cross correlation power 
spectrum for small £, leaving only the negative residual com- 
ponent of the TE power spectrum due to PGWs. Actually, 
this method can be treated a test of the (negative) contribu- 
tion to the TE correlation power spectrum due to PGWs on 
large scales using uncertainties in the measurements consis- 
tent with the total TE power spectrum. We use Monte Carlo 
simulations to analyze the probability of detecting PGWs 
using this method. 

By detection of PGWs we mean in this paper, the mea- 
surement of the parameter r, the ratio of the primordial ten- 
sor power spectrum, Pt(k), to the primordial scalar power 
spectrum, P s (k), taken at wavenumber, ko: 
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Pt(ko) 
Ps(k ) 



(1) 



where ko = 0.05 Mpc" 1 (see Smith et al. (2006)). For this 
paper we only consider the tensor-to-scalar ratio, r. All 
othe r parameters a re ass umed to be at their WMAP3 val- 
ues ( Spergel et al. ( 2007 )). The only other parameter which 
might affect the results is the tensor spectral index, nt, how- 
ever we assume, in this paper, that n t is very close to zero 
( Pciris et al, ( 200; )). The problem of n t < will be consid- 
ered in another paper. 

The plan of this paper is the following. In Section 
we describe the method for detection of PGWs based on 
measurements of the TE power spectrum. In Section ^, we 
describe the numerical Monte Carlo simulations we tested. 
In Sections ^, |j| and [], we introduce the statistical tests used 
to contstrain PGWs. In Section |^, we compare the three dif- 
ferent statistical tests used in this analysis. Section ra gives 



the r esults fo r the two toy experiments described in Polnarev 
et al. (jj007 ). The only uncertai nty in the first toy experi- 

1). In the second toy ex- 



ment is due to cosmic variance 
periment, along with cosmic variance, we take into account 
instrumental noise ([^) . We present results of Monte Carlo 
simulations for WMAP |0|) and Planck (Q. In Section [| 
we compare signal-to-noise ratio for BB and TE measure- 
ments. 



2 WIENER FILTERING OF THE TE CROSS 
CORRELATION POWER SPECTRUM 

Wiener filtering has been used often in the case of CMB 
data analysis. For example, it was used to combine multi- 
frequency data in order to remove foregro unds and extract 



the CMB signal from the observed data (Tegmark & Efs 



tathiou (1996); Bouchet et al. (1999)). Here we examine the 



use of the Wiener filter to subtract the PGW signal from the 
total TE correlation signal. This is done because the Wiener 
filter reduces the contribution of noise in a total signal by 
com parison with an estimation of the desired noiseless sig- 
nal (Vascghi ( 2006] )). In our case, the signal is the one due 
to PGWs only, and the signal contributed by density per- 
turbations is considered to be "noise" . 

The observed signal can be written as 
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where s and t refer to the contributions to the power spec- 
trum due to scalar and tensor perturbations respectively. 
The values aE.tm and aT,tm refer to the spherical harmonic 
coefficients of the temperature and polarization maps. In our 
application to TE correlation, we consider the Wiener filter, 
Wte.i' 
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The filtered signal, a' x tm (f° r X = T and E), is obtained 
from the measured signal, 
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In this paper, we assume the Wiener filter is perfect, in the 
sense that it leaves the signal due to PGWs only. We then 
get, for the filtered multipoles C k 
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In practice this is not true, because we are trying to 
determine Cf E , which is not known in advance. Neverthe- 
less, the assumption that the Wiener filter is perfect is good 
as a first approximation and illustrates the detectability of 
PGWs with the help of TE correlation measurements. 

The filtering can reduce the measured signal to the de- 
sired signal, but, since we are trying to remove the density 
perturbations and not the actual noise, we can not reduce 
the measurement uncertainties. These uncertainties in Cj E 
are then entirely determined by the noise in the original 
signal. 



From 
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020071 ), we showed that the TE 
power spectrum due to PGWs is negative on large scales, 
hence a test determining whether the Wiener filtered power 
spectrum is negative or not is a probe of PGWs. 

There are three different statistical tests we use to see 
if we can measure a negative TE power spectrum. The first 
test is a Monte Carlo simulation to determine signal-to-noise 
ratio, S/N (Section ^). The other two tests are standard non- 
parametric statistical test s: the sign test ( Section ^) and the 
Wilcoxon rank sum test ( |Wilcoxon| 01945| )) (Section |). 



3 MONTE CARLO SIMULATION 

For all of our tests, we calculate a random variable. If the 
data satisfies the hypothesis that r — 0, we can calculate 
the mean and uncertainty in the variables. If we make one 
realization of data, the random variable is determined from 
its distribution. Because we are not using any real obser- 
vational data, we must run a Monte Carlo simulation to 
reduct the risk of randomly getting a value for the variable 
taken from the outlying area of its distribution. To do this, 
the filtered multipoles, Cjf m , are randomly chosen from a 
gaussian distribution with mean Cjf and standard devia- 
tion ACT E where 
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(see, for example, Dodelson ( 2005 )), the variable f s k y refers 
to the fraction of the sky covered by observations and Nt 
is th e effective powe r spectrum of the instrumental noise 
(see Dodelson (2005 ) for details on how Nt is related to 
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Figure 1. This plot shows the combined distribution of all the 
Cff. 10, 000 values for every I in the range 2 $C I sC 53 were used 
in this histogram 



actual instrumental noise) . The under l ying power spectra 
are generated by CAMB[] ( [Lewis et al| ( 2000| )). 

Our determination of C e t is dependent on £. However, 
for two of our tests we ignore the value of £ in the calcula- 
tion of the random variable. We assume that the calculated 
random variable is gaussian. In order for this to work, the 
random variable must be calculated from gaussian variables. 
Fig. [j] shows the combined distribution of the multipoles for 
the "ideal" toy experiment (see Section pi). To do this, for 
every £ we determine 10, 000 different Cjj at values. We then 
combine all the different values, for every £, into one large 
distribution. Fig. |l] shows the errors on the multipoles are 
large enough so that for our statistical tests we can assume 
the multipoles are taken from a single distribution and not 
from a distribution that depends on £. 



4 MONTE CARLO S/N TEST 

For this test, the random variable we calculate, S/N, is de- 
fined as 



n TE 



(7) 



The reason why the sum in this equation is taken in 
the range 2 < £ < 53 is because only in this range 
sgn(Cf E (scalar)) = — sgn(Cj B (tensor)). In other words, 
if we include higher multipoles we confront with a danger of 
a false detection, because the total TE power spectrum is 
negative for £ > 53. 

The value of S/N is gaussian distributed because it is a 
sum of many modes of squares of gaussian distributed val- 
ues, Ce = a\ m . We approximate each Cj E as being gaussian 
distributed for the purpose of this paper. For each set of pa- 
rameters we run this simulation one million times to deter- 
mine the mean, (S/N), and standard deviation, &s/n- The 
mean of this distribution is determined by the preassumed 
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value of r, while the standard deviation is determined by 
parameters of the experiment and gives the confidence level 
of detection. We run such Monte Carlo simulations for dif- 
ferent values of r to determine in what range of r we can 
detect PGWs. When then using real observational data, we 
can compare the actual value of S/N with the results of 
Monte Carlo simulations to infer the likelihood, as function 
of r, which determines the probability that r / 0, or that 
PGWs exist at detectable levels. 



5 SIGN TEST 

The sign test is a test of compatability of observational data 
with the hypothesis that r = 0. If we do have r = 0, then 
Cjjm will be equally distributed around zero. Application 
of this test to the filtered data is very simple. In practice, 
all observational data are distributed between several bins 
and the averaging of the signal is produced in each bin sep- 
arately. Let Nf,i ns be the number of such bins. The sign test 
actually gives the probability that in JV_ bins the average 
is negative and in N+ = Nu-ns — N- it is positive, if r = 0. 
This probability, P, is given by the binomial distribution 
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The probability that the hypothesis r — is wrong is 
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The value X] ;=o ^ s ^ ne probability that we would get 
^ N + positive values given r = 0. This is the same as the 
probability of getting ^ N+ negative values given r = 0. 
Therefore our confidence that r 7^ is just 100% minus 
the sum of the probabilities describe above (the probability 
that the N+ is closer to the mean, Nu ns /2, if r = 0). This 
equation only makes sense if A+ < Nbms/2, since that is 
required for r > 0. If A+ > A(,, ns /2, that would imply 
r < 0, which is not physical. We would have to interpret 
the result as a random realization of r ^ 0, with the most 
likely result of r — 0. Therefore we would not be able to say 
r/not — with any confidence. 

Let us consider the following example: we put all mea- 
surements of Cj E into 11 bins and in three of them the 
average is positive. In this example, the probability that the 
hypothesis r = is wrong is equal to 89%. 

One possible drawback of this method is that it does not 
take into account any measure of the signal-to-noise ratio 
of individual measurements. As we show in Section (?], it is 
possible to have two completely different sets of data with 
the same probability of having r — 0. This test is also unable 
to make any prediction as to the value of r, only that it 
differs from zero. 



6 WILCOXON RANK SUM TEST 

This statistical test deals with two sets of data. The first set 
of data is taken from a real experiment which measures Cj E 
with some unknown r. The second set of data is generated 
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by Monte Carlo simulations (see Section ^]) with r = 0. 
The objective of the Wilcoxon rank sum test is t o give the 
proba bility that the hypothesis r = is wrong (Wilcoxon 
([1945|)). 



First, we choose some random variable U, whose proba- 
bility distribution is known if r = 0. For that, let us combine 
all data from first set with m multipoles and second set with 
7i2 multipoles into one large data set, which obviously con- 
tains n\ +?i2 multipoles. Then, we rank all multipoles in the 
large data set from 1 to n\ + 77,2 according to their amplitude 
(rank 1 for the smallest and rank n\ + 712 for the largest). 
Now, the variables Ri and R2 are defined as the sum of the 
ranks for the first original data set and the second original 
data set, correspondingly. Finally, the variable U, is 



U — min(J7i, U2), where 
Ui = Ri- ni (n, + l)/2, 



i = l,2 (10) 



If all multipoles of the first data set are larger than all mul- 
tipoles of the second data set, then U\ = n\n<x and U2 = 0. 
It is not difficult to show that Ui + U2 — niri2- If both sets 
of measurements have no evidence for PGWs, (Ui) = (^2). 
It is also simple to see that Ui + U2 — n\n?.. 

It is important to emphasize that the ranks of multi- 
poles are random variables because all multipoles themselves 
are random variables, hence Ui , U2 , and U are random vari- 
ables. If m + ri2 is large, the distribution of U can be ap- 
proximated as a gaussian with a known mean and standard 
deviation. In this approximation we have 
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In some cases, instead of U, the variables Ri or R2 are used. 
The reason U is used here is because mjj is symmetric in the 
data sets. If r = in both sets of data, then the distributions 
of Ui and U2 are the same, no matter what ni and 77,2 are. 
The distributions of Ri and R2 would be the same only if 
m =712. The probability that the first data set corresponds 
to r 7^ obtained from the test in which R\ or R2 is used is 
the same as if U is used. 

Since this test requires Monte Carlo simulations for the 
second set of data, we ran this test many times for many 
different data sets to get an accurate mean value for U. 

To reject the hypothesis r = means to detect PGWs. 
Using the Wilcoxon rank sum test the allowable value of r 
is determined, if instead of comparing with simulated data 
with r — 0, we compare with simulated data with r — ro 7^ 
0. In order to get a range of allowable values for r , we need to 
run multiple Monte Carlo simulations with multiple values 
for ro. This is where the assumption that the Cjf are from 
a random distribution that is independent of I is used (see 
Section ^). This implies that the ranks are random variables. 
If the errors on the Cff are small enough, then the ranks 
will be predetermined. Therefore, our assumption about the 
distribution of U will not be true and the test would have 
to be modified. Fortunately, this is not the case for even an 
experiment only limited by cosmic variance. 

To illustrate how this test works, let us consider the 
following example. Assume there are 4 multipoles in the first 



set of data and consider that r = 0.3 is the correct value. 
There are also 4 multipoles in the second set of data (which 
for sure corresponds to r = 0). All quantities below are 
expressed in /jK 2 . The value for the first data set are Ciq B = 



-0.005, Co, 



0.02, C 3 o = -0.015, and Cio = -0.01. 



The values for the second data set are Cj^ = 0.03, = 
0.003, C£f = -0.02, and Cf B = -0.003. A ranking of 
multipoles gives the ordering from lowest to highest, with 1 
referring to the first data set and 2 referring to the second 
data set, as 21112212. This results in i?i = 2 + 3 + 4+7 = 16, 
Ui = 16 - 10 = 6, and U 2 = 16 - 6 = 10. Therefore U = 
min{10,6} = 6. For m = 712 = 4, to reject the hypothesis 
that r = at 95% con fidence lev e l, Ui should be less than 
one (see, for example, Lehmann ( 1975 )). In this example, 
since Ui = 6 > 1, the first set of data cannot be considered 
as a detection of PGWs. 



7 COMPARISON OF TESTS 

The S/N test is greatly affected by outlying measurements. 
A measurement of one large negative multipole could falsely 
implay a detection. Both the sign test and the Wilcoxon 
rank sum test are not affected by individual outlying mea- 
surements. In the sign test, the value of individual measure- 
ments is irrelevant, because the test is sensitive only to the 
sign of individual measurements. The Wilcoxon rank sum 
test is affected by outliers, but considerably less than the 
S/N test. If the outlier is larger (or smaller) than every 
other multipole, its rank does not depend on its particular 
value. 

If we have two completely different sets of data, the 
main disadvantage of the sign test, as mentioned in Section 
II is that it could give the same result, while for the two other 
tests the chance to obtain the same value of r is negligible. 
For example, one set of data, consisting of 4 small negative 
multipoles and 4 large positive multipoles, gives the same 
result as another set of data, consisting of 4 large negative 
multipoles and 4 small positive multipoles. The S/N test 
gives two very different values of S/N for these two sets 
of data. We can also use the Wilcoxon rank sum test to 
compare these two sets of data. In this case U = 16 = \tuu, 
which corresponds to a confidence level of hypothesis that 
r = of less than 10%. 

With observational data, the sign test can be applied 
and does not require any Monte Carlo simulations (which 
could be considered as an advantage of this test). The S/N 
test requires Monte Carlo simulations, but only for the dis- 
tribution of the random variable S/N. The Wilcoxon rank 
sum test requires large Monte Carlo simulations and com- 
bines the data sets generated by these simulations with ob- 
servational data. In other words, Monte Carlo simulations 
are absolutely necessary after obtaining observational data, 
which may be considered a disadvantage of this test. Thus, 
each of the three tests has advantages and disadvantages, 
suggesting that the best way to work out observational data 
is to apply all these three tests. 
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8 DISCUSSION AND RESULTS 



The current best limit of r is r < 0.3 at 95% confidence pro- 
vided by WMAP in c ombination with previous experiments 
(spergel et al. (2007)). We need to see if this method can 
detect a value of r that is currently within the limit. We 
provide results for r = 0.3 in order to see how well these 
tests can detect an amount of PGWs that is very close to 
being ruled out by BB measurements. For all our tests we 
assumed that there is no foreground contamination. For the 
experiments that are not observing the full sky, correlations 
between multipoles must be taken into account. We bin to- 
gether the highly correlated multipoles so that the correla- 
tions between the bins are sufficiently small. 
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T he tw o toy experiments that were used in Polnarev 
( 2007 ) are also used here to c onstrain r. The t wo to y 
experiments are fully described in Polnarev et al. ( 2007 ), 
however we will also reproduce their description here. 

The first toy experiment is a full sky experiement. For 
this experiment, we take measurements over the full sky with 
no instrumental noise. The only uncertainty will be due to 
cosmic variance. This experiment represents the best limit 
to which the gravitational waves can be detected with the 
TE correlation. A space-based experiment with access to the 
full sky is the closest to this experiment. It is similar to what 
the Beyond Einstein inflation probe will be able to detect. 
This toy experiment will be hereafter referred to as the ideal 
experiment. 

The second toy experiment is a more realistic exper- 
iment. Measurements of the CMB are taken on 3% of the 
sky in one frequency (100 GHz) for 3 years. The noise in each 
detector of the 50 polarization sensitive bolometer pairs can 
be described by their noise equivalent temperature (NET) 
of 450 fiK^fs. The detector beam profiles are assumed to 
be gaussian and and it is described by their full width half 
maximum (FWHM) of 0.85°. 

This second toy experiment is similar to current ground- 
based experiments. Constraints from this experiment repre- 
sent those that can and will be obtained in the next several 
years using this method. This will be referred to as the re- 
alistic experiment. 

We also look at experiments with instrumental noise 
similar to the satellite experiments WMAP and Planck. The 
predicted errors for Planck are based on using the 100 GHz, 
143 GHz, and 217 GHz channel in the High Frequency In- 
strument (HFI). The numbers are gotten from the Planck 
science case, the "bluebook" Q The WMAP noise was ob- 
tained by using 3 years of the Q-band, V-band, and W-band 
detectors. 



8.1 Ideal Experiment 

A plot of the TE power spectrum due to PGWs for r = 0.3 
and nt = 0.0 for the ideal experiment is shown in Fig. [| 
The errors bars in Fig. ^ are calculated from the total TE 
power spectrum. There are no correlations between multi- 
poles, because we observe the full sky, but we plot the error 
bars binned in intervals of A£ =10 for simplicity in the plot. 
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Figure 2. The black line is the TE mode due to PGWs with 
r = 0.3. Red is the error bars for the ideal experiment calculated 
from the total TE power spectrum binned in intervals of A£ = 10. 





IdeaL experiment N + 



Realistic experiment N + 





Figure 3. This is a plot of the distribution of the number of 
positive multipoles for the Monte Carlo simulation for the ideal 
experiment (upper left), the realistic experiment (upper right), 
Planck (lower left), and WMAP (lower right). The dotted red 
line shows where N+ = i-/Vj, ins 



The Monte Carlo simulation gives an average of 19 mea- 
sured TE power spectrum multipoles greater than zero out 
of a total of 52 independent multipoles. If the null hypoth- 
esis was true, the sign test would indicate there is a 3.5% 
chance of measuring ^ 19 positive multipoles. This is equiv- 
alent to a ~ 1.8cr detection. A plot of the distribution of the 
number of positive multipoles is shown in the upper panel 
plot of Fig. ||. There is an 81% chance for the observed N+ 
to give a la detection of PGWs. 

The S/N test gives a mean value of S/N = —17.1 and 
standard deviation of 7.21. The upper left panel in Fig. [j] 
shows the distribution of the S/N values for the Monte Carlo 
simulation with r = 0.3. If r = 0.3 we would have a 0.8% 
probability of the measured S/N > 0. This negative value 
signifies that a non-zero tensor-to-scalar ratio produced an 
anti-correlation. We can assume that the standard deviation 
would be the same if the mean of S/N was (equivalent to 
r — 0.0), because it is equivalent to adding a constant value 
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Figure 4. The S/N statistic distribution for the ideal experiment 
(upper left), realistic experiment (upper right), Planck (lower 
left), and WMAP (lower right). The dotted red line shows where 
S/N = 0. 



Figure 6. This is the plot of the distribution of U for the ideal ex- 
periment (upper left), realistic experiment (upper right), Planck 
(lower left), and WMAP (lower right). The red dotted line is the 
value for mjj and the light blue dashed lines enclose the ltr region 
for U assuming the hypothesis that r = 



10 




-30 E i 

0.0 0.1 0.2 0.3 0.4 0.5 



r 

Figure 5. This is a plot of (S/N) and crg/ N as a function of r 
for the ideal experiment. The black line is {S/N) and the red line 
is cr s/N . 



to every measured value (and hence adding a constant to 
S/N which would not change the error). Therefore, if r — 0, 
the probability of getting S/N < —17 .4 is 0.8%, and hence 
we have a 99% chance that A plot of (S/N) and as/N 

as a function of r is shown in Fig. H. As can be seen from 
the plot, we can predict a value of r for any value of S/N . 
The value of (Js/n is a relatively constant function of r and 
so our prediction about the distribution of S/N for different 
value of r is a good approximation to the true distribution. 

The Wilcoxon rank sum test gives Uavg — mu = 
— 1.23<7Lr • The variable Uavg is the mean value for U in the 
Monte Carlo simulations described earlier. The values mu 
and au are given in Section ^| The distribution of U for 
the Monte Carlo simulations with r = 0.3 is shown in Fig. 
^[ The standard deviation of the distribution of measured 
U is the same as the standard deviation of the distribu- 
tion of U assuming the hypothesis that r = 0. The only 
difference between the distributions is that mu is shifted 
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Figure 7. This is the plot of the signal-to-noise ratio (number 
of as) for different values of r for the three different tests. The 
black line is the S/N test, the dashed dark blue line is the sign 
test, and the dotted-dashed light blue line is the Wilcoxon rank 
sum test 



by a constant value. Therefore, there is a 22% chance that 
U — mu < —2ou- There is also a 40% chance that we mea- 
sure U — mu < —lau, and are not even able to make a la 
detection of PGWs. 

A comparison of the three tests is shown in Fig. |^. This 
is obtained by simulated with with several values of r and 
then interpolating between them. A 2a detection is obtained 
for r = 0.26 (S/N test), r = 0.3 (sign test), and r = 0.5 
(Wilcoxon rank sum test), highlighting its intended use as 
a monitor of a false positive detection for large r. 

8.2 Realistic Ground Based Experiment 

A plot of the TE power spectrum due to PGWs with r = 0.3 
is shown in Fig. H. The error bars are calculated from the 
total TE power spectrum. Observations of an incomplete 
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Figure 8. This figure is the same as Fig. |i| except for the realistic 
experiment 



sky require the multipoles to be binned in sizes of A£ = 10. 
This experiment has much larger error bars than the ideal 
experiment and it will not be able to detect low values of r 
with the TE cross correlation only. 

The results for the Wiener filtering method were much 
worse than those for the ideal experiment for r — 0.3. Since 
this experiment observes a small portion of the sky, the mul- 
tipoles are correlated and we must bin together to get rea- 
sonably uncorrelated measurements. For this experiment, we 
only have 7 to 8 uncorrelated multipoles, instead of 52 uncor- 
related multipoles in the case where the full sky is observed. 
Getting 7 out of 8 negative multipoles is a 3% probability if 
there are no PGWs. For the Monte Carlo simulations of the 
realistic experiment, on average, half of measured multipoles 
are positive and half are negative. A plot of the distribution 
of the number of positive multipoles is shown in the upper 
right panel of Fig. ^. In this case, we cannot distinguish 
r — 0.3 from r — 0.0 with any significance. 

The S/N test gives an average value of S/N = —0.95 
with standard deviation of 2.64. For the realistic toy exper- 
iment, the distribution of S/N for r = 0.3 is shown in the 
upper right panel of Fig. ^. In order to obtain 68% confidence 
detection of PGWs, we must use r ~ 0.7. In this sense the 
TE test provides monitoring and insurance against false pos- 
itive detection with r > 0.7, which could arise, for example, 
if foregrounds or other systematic effects arer improperly 
removed. 

The last statistical test, the Wilcoxon rank sum test, 
gives Uavg — rnu = — 0.20au- The distribution of U for 
r = 0.3 is shown in the upper right panel of Fig. ^| This 
gives the weakest result in terms of the three tests for the 
Wiener filtered data. The realistic experiment will not be 
able to constrain r < 0.3 using the TE cross correlation 
power spectrum. Its limit is closer to r < 0.7 — 0.9 at only 
68% confidence depending on the test used. For a higher 
confidence in a detection of PGWs, the value of r would 
need to be much higher. Since the observed distribution of 
U corresponds almost exactly to the simulated distribution 
of U under the assumption that r = 0, therefore we have a 
16% chance of measuring U — mu < —\ajj. 



The results of the Wiener filtering showed that the WMAP 
cannot make a detection of gravitational waves using the TE 
cross correlation power spectrum alone. As with the two toy 
experiments, the result of the scalar and tensor separation 
was similar. The Monte Carlo simulation gave on average 
gave 13 positive multipoles out of a total of 26 uncorrelated 
multipoles. We would get the same result if the input data 
had r = 0.0 so we cannot detect PGWs with WMAP using 
only the TE power spectrum. A plot of the distribution of 
the number ofpositive multipoles is shown in the lower right 
panel of Fig. j3|. As can be seen, this distribution of 7V+ for 
WMAP noise and r = 0.3 is simply the distribution for 
r = 0. 

For WMAP, the S/N test gives the value of S/N = 
—0.02 with a standard deviation of 5.09. The distribution is 
shown in the lower right panel of Fig. ^| The distribution 
is centered around S/N — so there is no chance of using 
this test to detect PGWs in WMAP's TE power spectrum. 
The probability of getting a la or 2a detection is the same 
probability that we would randomly get a detection if there 
are no PGWs. 

The rank sum test gives a value of Uavg — mu = 
—0.004(7(7, which is implies no ability to distinguish 
WMAP's observed TE data from a data set with no PGWs. 
A plot of the distribution of U for WMAP error bars is shown 
in the lower right panel of Fig. [| We reach the same conclu- 
sion for WMAP noise as for the realistic experiment. There 
is only a 16% chance that we can measure U — mu < — lau 
and make a la detection of r = 0.3 

The published WMAP results show an anti-correlation 
of TE power spectrum at large scales. Un fortunately this 
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i s not a detection of PGWs as theorized in 
( |2006| ). The contribution to the TE power spectrum due 
to PGWs only changes sign once for I < 90. If a claimed 
evidence for gravitational waves is to be believed, then the 
TE power spectrum would have to change sign three times 
for I < 60. In fact, other than the two anticorrelations at 
low I, the rest of the multipoles, up to I = 53, are consistent 
with r = 0. None of the described tests applied to the current 
WMAP data will give any detection of PGWs. 



8.4 Planck 

The sign test gives on average 10 positive measurements of 
the TE power spectrum out of a total of 26 uncorrelated 
multipoles. There is a 16% chance of getting ^ 10 positive 
multipoles if r — 0. A plot of the distribution of the number 
of positive multipoles for Planck is shown in the lower left 
panel of Fig. || There is a 50% chance that we will measure 
N+ < 10 and hence have a la detection of r = 0.3. 

The S/N test gives a value of S/N = -6.24 with a 
standard deviation of 5.09. There is only a 10% chance that 
the S/N test results in a value of S/N larger than zero, if r = 
0.3, and a 10% chance getting S/N < -3.12 if r = 0. This 
is close to a 90% probability of detection. The distribution 
of the S/N variable is shown in lower left panel of Fig. 

Again, the rank sum test gives the lowest confidence 
result with a value of Uavg — rnu = —0.66au- A plot of the 
distribution of U is shown in the lower left panel of Fig. ^[ 
There is a 37% probability that we will measure U — mu < 
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-lau and a 9% probability that we measure U 
-2au for Planck. 



mu < 



9 COMPARISON OF MEASUREMENTS OF 
TE POWER SPECTRUM WITH BB POWER 
SPECTRUM 

CMB polarization can be separated into two distinct compo- 
nents: E-mode (grad) polarization and B-mode (curl) polar- 
ization. PGWs generate B-mode polar ization i n con t rast to 



density perturb ations (see for example Seljak ( 1997 ) ; Bcljak 



fc Zaldarriaga (1997); Kamionkowski & Kosowsky (1995)) 
Therefore, most CMB polarization experiments searching 
for evidence o f PGWs focu s on m e asuring the BB powe r 
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spectrum (s ee Taylor et al. (2004) 



Bowden et al 



(2004); 



(2006)). For different aspects of CMB polariza- 



t ion g en erated by PGW s s ee, for example, Basko fc Polnarev 
([1980 ) ; |Polnarev| (|l985h ; JCrittenden et al] (|1993T 



k Kosowsky] ( p98[ ) 
( [200(:S| ) 



Seljakl (|1997|) ; IScljak fc Zaldamaga| ( |1 99 
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Frewin 



et al.| f[1994D; |Coles et aL| (|1995|); fcamionkowski et al | (|1997[) 



Kamionkowski 
]Kelrtmg'*e't*ld 



We have shown that the TE cross correlation power 



spectrum off e rs an other method of detecting PGWs (Grit 
(1995)). The TE power spectrum is two orders 
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of magnitude larger than the BB power spectrum and it was 
originally suggested that it might be easier to detect gravi- 
tational waves in the TE power s p ectrum than using the B B 
power spectrum (Baskaran et al. ( 2006 ); Grishchuk| (| 2007 ) ) . 



However, as we have shown in Polnarev et al. (20071), that, 
when we use zero multipole method, uncertainties in TE 
measurements exceed uncertainties in BB measurements in 
such a way that the signal-to-noise ratio for BB measure- 
ments is better than for TE measurements. We show below 
that the same is true for the separation of scalars and ten- 
sors. 

However, an advantage of TE measurements for ground 
based experiments, which observe only a small fraction of 
the sky, is related to the fact that the main spurious effects in 
the BB power spectrum are cuase d by E/B mixing. This w ill 
limit the r that can be detected (Challinor & Chon (2005)). 
The E-modes are practically unaffected by E/B mixing, so, 
in contrast to the BB measurements, the TE power spec- 
trum should be nearly the same for both full and partial sky 
measurements. 

Along with this, the methods based on the TE cross 
correlation can be considered as very useful auxiliary mea- 
surements of PGWs, because systematic effects in TE mea- 
surements are independent from those in BB measurements. 
For example, T/B leakage or even E/B leakage could swamp 
a detection of BB, w hereas T/E lea k age w ould be small and 



well controlled (see Shimon et al 



( |2007| )). These BB sys- 
tematics could falsely imply a large r, but measurements 
of the TE power spectrum provide insurance against such 
a spurious detection. Additionally, galactic foreground con- 
tamination will affect BB and TE in different ways, which 
enable us to perform powerful cross-checking and subtrac- 
tion of foregrounds in BB measurements. 

For the zero multipole method, we explained why the 
signal-to-noise ratio for BB measurements is better than for 
TE measurements. Below we give simple summarizing argu- 



ments why the same is true for the Wiener filtering of the 
TE power spectrum. 

If Nf <C Cf B , the signal-to-noise ratio for the BB power 
spectrum is 



(S/N) BB = 



r <BB 

AC BB 



= 7 



r <BB 



C BB + Ni 



where 



7 : 



(2e + i)f sky 



(13) 



(14) 



If N e > Cf B then we will not be able to detect PGWs and a 
comparison with the TE power spectrum is not worthwhile. 

If N t < Cf E and r < 1, for the TE power spectrum, 
the signal-to-noise ratio is 



(S/Nh 



ACf E 
V2-) 



[{CJ E Y + {CJ T + N e /2){C EE + Nt)] 1 ' 2 



V2-Y: 
V2-y 



C 



[(Cf E ) 2 + Cf T C EE ] 1/2 



a + f3r 



(15) 



where a and (5 are 



P = 
where 



DJ E 

nr<TB n TE i r^TT n EE . s~iTT r\EE 

2DT E a ' 



r,XY r ,XY I 



(16) 



(17) 



One can see that a and j3 are on the order of unity. There- 
fore, the signal-to-noise ratio is approximated as 



{S/N) TE = s/2-y 



a + (3r 



V2 7 - 



(18) 



In other words ifr<a//3~l, BB measurements have the 
obvious advantage in comparison with the Wiener filtering 
of the TE power spectrum. Indeed if r < 1, (S/N)bb ~ 7, 
while (S/N)te ~ jr < 7. This is because in BB measure- 
ments, applying proper data analysis, we can entirely elim- 
inate contributions of scalar perturbations to CMB polar- 
ization signal as well as to the uncertainties. For the perfect 
Wiener filtering of the TE power spectrum, we can eliminate 
the contribution of scalar perturbations to the signal only, 
but cannot eliminate their contribution to the uncertainties. 



10 CONCLUSION 

The method described here is one in which we filter out 
the signal due to density perturbations, leaving only the 
contribution to the TE power spectrum due to PGWs. We 
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then test the resulting TE power spectrum to see if it is 
negative. Three different statistical tests were used to see if 
there was a significant detection of PGWs. The S/N test 
can give a value for r using a comparison with Monte Carlo 
simulations, while the Wilcoxon rank sum test can only give 
an allowable range for r. The sign test will only tell us if 
r / 0. 

From Polnarev et al. (2007), we saw that we could de- 
tect r = 0.3 to 3a with a measure of to, the position where 
the TE power spectrum first changes sign, for the ideal ex- 
periment. Using the method discussed in this paper, we see 
that we are unable to make this significant of a detection. 
The best result was for the S/N test which would give a 
2.3a detection of r = 0.3. To detect PGWs on the level of 
3a, the tensor-to-scalae ratio r should be r ^ 0.4. The sign 
test would give 2a detection for r — 0.3 and a 3a detection 
for r — 0.45. The Wilcoxon ranked sum test gives only a 
1.2a detection for r = 0.3 and a 3a detection for r — 0.7. 
Similar results were gotten for the other three experiments 
tested. Thus in the sense of potential to detect PGWs, the 
zero multipole method is the best, next best is the S/N test, 
then the sign test, and the worst is the Wilcoxon ranked sum 
test. 

Baskaran et al. ( ^006| ) present illustrative examples in 
which high r is consistent with measured TT, EE, and TE 
correlations. The value of r is so high in these examples that 
if PGWs with such r really existed, current BB experiments 
would already detect PGWs. All models predict that the TE 
cross correlation power spectrum change sign only once for 
I < 100. The fact WMAP cannot exclude several multipoles 
with Cf B > in between multipoles of Cj E < means that 
the TE cross correlation power spectrum either changes sign 
several times for I < 100 or there is some instrumental noise 
which causes some anticorrelation measurements. Using in- 
strumental noise consistent with WMAP, our Monte Carlo 
simulations give A^o ~ 16 and £q > 40, which means that 
there is no evidence of PGWs in the TE correlation power 
spectrum. 
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